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Abstract. We discuss the pairing gap, a measure for nuclear pairing correlations, in chains of spherical, 
semi-magic nuclei in the framework of self-consistent nuclear mean-field models. The equations for the 
conventional BCS model and the approximate projection-before-variation Lipkin-Nogami method are for- 
mulated in terms of local density functionals for the effective interaction. We calculate the Lipkin-Nogami 
corrections of both the mean-field energy and the pairing energy. Various definitions of the pairing gap are 
discussed as three-point, four-point and five-point mass-difference formulae, averaged matrix elements of 
the pairing potential, and single-quasiparticle energies. Experimental values for the pairing gap are com- 
pared with calculations employing both a delta pairing force and a density-dependent delta interaction in 
the BCS and Lipkin-Nogami model. Odd-mass nuclei are calculated in the spherical blocking approxima- 
tion which neglects part of the the core polarization in the odd nucleus. We find that the five-point mass 
difference formula gives a very robust description of the odd-even staggering, other approximations for the 
gap may differ from that up to 30% for certain nuclei. 

PACS. 21.60.Jz - 21.30.Fe - 21.60.-n 



1 Introduction 

Pair correlations, which play a crucial role in supercon- 
ducting solids [GJ, also constitute an important comple- 
ment of nuclear shell structure Most often, pairing is 
treated in the so-called BCS approximation which 
was introduced in the original paper of Bardeen, Cooper 
and Schrieffer Q . The nuclear applications are particular 
in two respects: first, nuclei are finite, in fact rather small, 
objects, and second, we do not yet have a sufficiently re- 
liable microscopic nuclear many-body theory from which 
we could deduce the appropriate pairing interaction and 
its strength. The second problem causes two further ques- 
tions. First, one needs to develop a reliable and manage- 
able form for the pairing energy functional, and second, 
one has to determine an appropriate pairing strength for 
a given functional. 

There exist various prescriptions for the pairing en- 
ergy functional. Schematic pairing forces basically consist 
of defining a small band of pairing-active states and pa- 
rameterizing one typical pairing matrix element in depen- 
dence on the system size, i.e. neutron and proton num- 
ber. They are convenient and successful in many respects. 
But they are plagued by serious disadvantages: the cou- 
pling to continuum states is much exaggerated and the 
parameterization in terms of system size becomes ques- 
tionable for deformed systems along the fission path. To 



avoid these problems local two-body pairing forces are in- 
creasingly used which is particularly satisfying in connec- 
tion with self-consistent mean field calculations 
There are even some mean-field models like the Gogny 
forces [|l0]|ll| or the particular Skyrme force SkP |l2| which 
aim at a simultaneous description of the particle-particle 
and particle-hole channels of the effective interaction us- 
ing the same force. That is not a necessary condition. 
The simple local pairing energy functionals also provide a 
very good description of pairing properties throughout the 
whole chart of isotopes using only two universal strength 
parameters, one for protons and one for neutrons. This 
makes these forms for the pairing energy functional more 
reliable for calculations of deformed nuclei and of nuclei far 
off the valley of stability than the widely used schematic 
pairing forces. In the following we will concentrate on the 
class of local pairing interactions. 

The small particle number of nuclei often interferes 
with the fact that the BCS ground state mixes particle 
numbers with a relative spread of order 1 / \/~N. In princi- 
ple, one has to perform a projection of the BCS state be- 
fore variation, but this exact particle number projection 
can be very cumbersome, see, e.g., |l^,Q. The Lipkin- 
Nogami (LN) scheme offers a reliable approximate projec- 
tion method jl5| 16 , l?],[l8| . It has the technical advantage 
that it is formulated completely in terms of BCS expec- 
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tation values which makes its numerical implementation see e.g. [12 



very simple, for detailed discussion see |17| , [18| . [19[|2(| . We 
will discuss the LN scheme side by side with the BCS 
treatment. 

There remains as the last and crucial problem the de- 
termination of an appropriate pairing strength. Insuffi- 
cient microscopic information requires that one recurs to 
a phenomenological assessment. This line of development 
has been followed with increasing accuracy since the intro- 
duction of pairing in nuclei, for a comprehensive compila- 
tion of pairing with schematic forces see . Also the local 
two-body pairing interactions leave the overall strength 
as a free parameter to be fixed phenomenologically. This 
has to be done with respect to an observable sensitive to 
pairing correlations. Such an observable would ideally be 
provided by the pairing gap which, however, is not di- 
rectly accessible in experiments. The observable quantity 
with probably closest relation to the pairing gap is the 
pronounced odd-even mass staggering of nuclei, which is 
usually used for fitting the pairing strength. Even here, 
though, there remains a choice of several recipes on how 
to extract the pairing gap or odd-even staggering, respec- 
tively, from a combination of neighboring mass values. 
From the theoretical side, of course, one has more direct 
access to a gap. But even there ambiguities emerge and 
one is left with several possibilities to define measures for 
the pairing correlations. 

It is the aim of this paper to compare and discuss 
the various definitions of a pairing gap in models of lo- 
cal two-body pairing interactions in connection with self- 
consistent mean-field approaches, both at the level of the 
conventional BCS approach as well as of the LN method. 
Furthermore, we will compare two variants of the local 
zero-range pairing forces, a delta force and a density- 
dependent form |9[[l0|. Last but not least, the discussion 
will extend to exotic nuclei where differences between the 
various definitions for the pairing gap and options for the 
pairing method become particularly obvious. 

The paper is outlined as follows: In Section[2]we present 
the basic equations of the Lipkin-Nogami model employ- 
ing local interactions in the framework of the Skyrme- 
Hartree-Fock model needed for our discussion. In Sec- 
tion H various approximations for the pairing gap are dis- 
cussed, in Section |] the results are presented, Section || 
summarizes our findings. 



2 The Theoretical Framework 



and such nuclei are not considered here. All 
conclusions about pairing gaps drawn in this paper can 
thus safely be derived in the BCS approximation. 



2.1 The Mean Field 

Presently the most widely used self-consistent mean-field 
models are the non-relativistic Hartree-Fock approach with 
either the Skyrme (SHF) fH or the Gogny force 0, 
and the relativistic mean-field model J23|. They all pro- 
vide a well-adjusted effective energy functional for nuclear 
mean-field calculations. We choose the SHF model for the 
present investigation. The description starts from an en- 
ergy functional 



f — f 



(1) 



whose mean-field part £ m f is formulated in terms of the 
local distributions of density p, kinetic density r, and spin- 
orbit current J. At this point it is not necessary to unfold 
all details of this rather elaborate functional, we abbrevi- 
ate the dependence with the most general case, the full 
one-body density matrix 



p= P (x,x') = (v; t (x')vi(x)) 



(2) 



from which all loc al d ensities and currents can be de- 
rived, see appendix |A.2| for details, ffi (x) creates a particle 
with spin projection a/2 at the space point r. Through- 
out this paper (■ ■ • ) denotes BCS expectation values. For 
the Skyrme energy functional we choose the rather recent 
parameterization SkI4 p4| . The pairing energy functional 
depends addition ally on the local pair density \ as intro- 



duced in Section |2.2| . Variation of the energy functional 
£ with respect to the single-particle wave functions 
yields the mean-field equations 



h<t>k = £fc 4>k with h = — 



S£ 
Sp 



(3) 



We investigate the pairing gap in the framework of self- 
consistent mean-field models. Pairing correlations are treated 
on the HF+BCS level where the equations of motion are 
derived by independent variation with respect to single- 
particle wave functions and occupation amplitudes. This 
is a widely used approximation to the more involved HFB 
approach where wave functions and occupation amplitudes 
are varied simultaneously [^j. The BCS approximation is 
applicable for all well-bound nuclei, i.e. for most of the 
nuclei discussed throughout this paper. It becomes criti- 
cal only for nuclei close to the neutron or proton drip line, 



where we have neglected the contributions from the vari- 
ation of the pairing density x i n the energy functional to 
the equations-of-motion of the single-particle states <j>k- 
This constitutes the BCS approximation to pairing (see 
[P5[ for the discussion of the full HFB equations in the 
representation in natural orbitals that is used here). The 
mean-field equations are solved on a grid in coordinate 
space with the damped gradient iteration method and a 
Fourier representation of the derivatives. The numerical 
techniques are summarized in |26| . 



2.2 The Pairing Energy Functional 

We parameterize the effective pairing interaction in terms 
of a local pairing energy functional of the form 



E 



d 3 r X *(r) Xq (r) G,(r) 



(4) 
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which allows for a spatial modulation of the strength G(r). 
x(r) is the local part of the pair density matrix 



X q (r,a;r,a) 



E 

<T=± 



= -2 



E 

k£S2 q 



Uk Vk |0fc(r)] 5 



a {i) q {v,-a)i) q {v,a)) 



(5) 



the Fermi energy. AE q is fixed from the condition that 
the sum of the cutoff weights includes approximately one 
additional shell of single-particle states above the Fermi 
surface 



E h = N q + i.mN t 



2/3 

q 



(9) 



ken„ 



with q G {p, n}. The 4> k are the single-particle wave func- 23 The Lipkin-Nogami Equations 

tions and Vk, Uk = — v\ the pairing amplitudes. We 



restrict ourselves to stationary states of time-reversal in- 
variant systems and pairing between like particles only. 
This sorts the single-particle states into conjugate pairs 
k «-> k = — fc with the same spatial properties but opposite 
projection of the total angular momentum and allows to 
restrict the summation in the pair density to k > 0. The 
time-reversal symmetry also renders the pair density real, 
i.e. x*(r) = x(r). 

Two models for the spatial modulation of the pairing 
strength are considered here 



The LN scheme serves as an approximation to particle 
number projected BCS. It can be derived by a mo mentum 
expansion of the projected BCS equations 15 111" 



At the end this boils down to adding overlaps with the 
variance of the particle number (AN q ) 2 at various places. 

In most cases, the LN method is used with a simple 
schematic pairing interaction in the framework of macro- 
scopic-microscopic models and self-consistent models for 



eneraysurfaces Plflp^,|33|,p 


35 






43 


HggTOnly re- 



G g (r) 




P(r) 
Po 



DF, 
DDDI. 



as well as high-spin states 

cently, theLN scheme was employ ed for a local delta pair 
ing force 



(6) 



Usually only the 
culated; but in self- 



The simpler case (DF) can be deduced from a delta force 
for the pairing interaction Vp a ir(r, r') = Vo, q S(r - r'), 

while the other (DDDI) corresponds to a density -dependent 



delta interaction [p|27 



, 2q| . The additional parameters of 
the DDDI force are set here to 7 = 1 and po — 0.16 fm~ 3 
(i.e. the saturation density of symmetric nuclear matter). 
More general choices are conceivable Jl^j2^] but very hard 
to adjust phenomenologically. Thus we keep to the sim- 
plest choice above. Note that a separate pairing strength 
Vo, p or Vb,n is associated to each nucleon sort. This ex- 
plicit breaking of the isospin symmetry in the pairing en- 
ergy functional is standard in nearly all pairing forces and 
schematic models, see e.g. |2l| , |30| . 

Although the pairing matrix elements deduced from 
the pairing functional (0) suppress the contribution from 
unbound states located outside the nucleus considerably 
(as compared to the schematic pairing force), the implicit 
zero-range nature of the pairing force still tends to overes- 
timate the coupling to continuum states. This defect can 
be cured to some extent using finite-range forces like the 
Gogny force [jlO 1 1 which are, however, cumbersome to 
handle. We prefer to simulate the effect of finite range by 
introducing smooth energy-dependent cutoff weights 18] 



J and the Gogny force 1 47 46 
correction of the pairing energy is caL 
consistent models the contribution of the mean field to the 
total binding energy is calculated from the BCS state as 
well, so that the correction of the pairing energy has to 
be complemented by a correction of the mean-field energy 
as considered in p8| , p7p^ |. In this paper, we present and 
employ the LN equations in the context of self-consistent 
mean-field models and for the case of local pairing energy 
functionals. As pairing gaps are the theme of this paper, 
particular emphasis is laid on the properties of the LN 
scheme relevant for the discussion of pairing gaps. 
The LN equations are derived by variation of 



K = £ - E (W^)+A 2 , 9 (7V 9 2 ) 

q£p,n 



(10) 



Variation of ( |i0| ) with respect to the occupation ampli- 
tudes Vk leads to 



1 - 



VWk-\) 2 + f^l 



(11) 



This is the standard expression for the occupation num- 
ber v 2 . in the BCS model ||||[3), nere containing a state- 
dependent single-particle gap A k times the cutoff factor 
fk and a generalized Fermi energy 



fk — 



1 



\ q = A l!9 + 4A 2i9 (7V 9 + l) 



(12) 



1 + cxp[(e fe - X q - AE q )/fi q ] 
in the evaluation of the local pair density 

f k u k v k |<fe(r)| 2 



Xq( 



E 

f!>0 



(7) 



(8) 



which is determined from a constraint on particle number. 
The quantity e' k is a renormalized single-particle energy 



4 = e -k + 4 \ 2 , q v k 



(13) 



In case of time-reversal invariance, the state-dependent 
single-particle gaps are given by 



The cutoff parameters AE q and p q = AE q /10 are chosen 
self-adjusting to the actual level density in the vicinity of 



A k 



d'r <f>l(r)A q (r)4, k (r) 



(14) 
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i.e. they are matrix elements of the local pair potential 
A?( r ) = |fe = ^X q (r) G„(t) . (15) 



Note that X 2 is not a Lagrange parameter |l^ , |33| . It is 
held fixed during the variation and is determined after 
variation from the additional condition |Ll7fl 



^2,q — 



((H nl{ + H pail )(AN 2 , q ) 2 ) 



(16) 



For simplicity of the presentation, Eq. ( |l6| ) is written for 
the case of an underlying many-body Hamiltonian H . The 
discussion of the more general case of an energy functional 
used in this paper is presented in Appendix A.l . N 2 ,q is the 
part of the particle-number operator that projects onto 
two-quasiparticle states 



fee O t 



u k v k {a\a\ + ai& k ) 



(17) 



while AN 2 . q = N 2 , q - (N 2tq ). The numerator of Eq. (|16|) 
contains, besides the familiar contribution from the pair- 
ing functional, an additional one from the linear response 
of the mean-field to the particle-number projection, see 
[ p8| . The total binding energy after approximate particle- 
number projection is given by 



E 



LX 



= £- E A ^ 

<?Gp,n 



(18) 



For arbitrary one-body operators the LN expectation value 
can be calculated introducing effective LN occupation num- 
bers and local densities, see [ [33||38| . 

Thus far the presentation applies to the more involved 
LN method. The BCS approximation is recovered simply 
by setting X 2 ^ q — in the above equations. 

2.4 Blocking 

The evaluation of the odd-even staggering involves also 
nuclei with odd mass number where one pair of conju- 
gate states has to be blocked, i.e. taken out of the pair- 
ing scheme. One of the blocked states has the occupation 
Ufc = 1, the other vj, — 0. In the standard textbook ap- 
proach the blocked many-body state is constructed non- 
self-consistently from the BCS ground state as a one- 
quasiparticle excitation, see e.g. [0Jq,p] and Section 3.3. 



The generalization of this approach for energy functionals 
is outlined in Appendix [b|. 

In a self-consistent approach to the blocked many- 
body state the single-particle wave functions and occupa- 
tion amplitudes have to be determined from a variational 
principle. Blocking a pair of states changes the density 
matrix, Eq. (0) and with that the single-particle Hamil- 
tonian, Eq. (ph. Besides a rearrangement of the single- 
particle wave functions the unpaired nucleon causes a po- 
larization of the core by breaking rotational and time- 
reversal invariance in the intrinsic frame (we assume spher- 
ical BCS ground states only), see e.g. [|49| . This requires a 



deformed calculation of the odd-mass nucleus considering 
also time-odd contributions to the single-particle Hamil- 
tonian (||) as discussed in [p9| , ^o|| . 

The change in binding energy due to the core polariza- 
tion depends on the properties of the time-odd spin and 
spin-isospin channels of the effective interaction which are 
not yet well adjusted for current mean-field models, see 
e.g. |51[] and references therein. Calculations with the cur- 
rently available models suggest that the polarization is 
non-ne gligibl e for the description of the odd-even stagger- 
ing p9| , p2|5^ , |5^ , p5| , but effective interactions with prop- 
erly adjusted spin and spin-isospin channels are needed 
before the effect can be treated quantitatively. In view of 
these uncertainties we restrict ourselves here to the sim- 
pler spherical blocking approximation, where one replaces 
the blocked single-particle state by an average over the 
degenerate states in its j shell, restoring rotational and 
time-reversal invariance of the many-body system in the 
intrinsic frame. In practice this means that the weight of 
the blocked j shell is given by (2j — l)uv when calculating 
the pair density and (2j — l)v 2 + 1 when calculating the 
local densities and currents. All other states enter with 
their full degeneracy 2j + 1. This approximation includes 
the large part of the rearrangement effects from monopolc 
polarization, but omits the polarization effects from mul- 
tipole deformations and time-odd currents. 

Owing to the rearrangement effects blocking of the 
single-particle state with smallest quasiparticle energy (as 
defined in Sect. 3.3) does not necessarily lead to the largest 
possible total binding energy. Therefore one has to per- 
form calculations for a number of blocked single-particle 
states around the Fermi energy and search for the config- 
uration giving the largest total binding energy. 



3 The Pairing Gap 

3.1 Nuclear Masses and Odd-Even Staggering 

The key feature of pairing correlations is the occurrence of 
an energy gap in the excitation spectrum. This gap man- 
ifests itself in two different kinds of energetic observables: 
First, there is a gap in the quasiparticle excitation spectra 
of even-even nuclei, which does not appear in the spec- 
tra of odd-mass number or odd-odd nuclei, and second, 
there occurs a shift between the interpolating curves of 
the ground-state binding energies of even-even as com- 
pared to odd-mass nuclei, which is called the odd-even 
mass staggering. Usually, the second phenomenon is ex- 
ploited to define the experimental pairing gaps assuming 

M 



E eve n-even(Z , N) — E (Z,N) , 

E odd z(Z, N) = E {Z, N) + A P (Z, N) 
E odd N (Z, N) = E (Z, N) + A n (Z, N) 



(19) 



In odd-odd nuclei there is additionally the residual inter- 
action between the unpaired proton and neutron, but this 
case will not be considered in the present discussion. In a 
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self-consistent mean-field approach Eq is the (negative) 
energy of the fully paired many-particle wave function, 
i.e. the BCS ground state, while A q is the energy lost due 
to the blocking of a pair of conjugate states by the odd nu- 
cleon. The gap introduced with (|l^) has to be interpreted 
carefully The gap as defined through the separation ( fl9| ) 
contains more than pure pairing correlations. It includes 
unavoidably all p olar ization effects from the mean field as 
outlined in Sect. 2.4 . Despite of these uncertainties ( |l9|) 
provides the starting point for the definition of experi- 
mentally accessible pairing gaps which will be discussed 
in Sect. ^2. The gap defined with ( |i~^ ) serves than as a 
point of reference, as it can be calculated directly within 
the mean-field model 

A<- h \Z,N) :=E hlock {Z 1 N)-E {Z 1 N) (20) 

as the difference in binding energy between the one blocked 
state of an odd-mass number nucleus with the largest 
binding energy and its fully paired (fictitious) BCS vac- 
uum. Both calculations have, of course, to be performed 
self-consistently. We will call A^ b ' the "blocking gap" in 
the following. 



3.2 Finite-Difference Formulae 



The gap from Eq. (20) is a purely theoretical construct. 
The problem is that Eq is not measurable for odd-mass 
nuclei. We need a definition which is also experimentally 
accessible. It should fulfill two requirements which are use- 
ful for the phenomenological adjustment of the pairing 
strength: first, it should be easy to calculate both theoret- 
ically and from experimental data, and second, it should 
be influenced as little as possible by the properties of the 
underlying mean field in order to decouple mean-field and 
pairing properties. The odd-even staggering is related to 
differences of nuclear masses and as such easy to measure 
as well as to compute, As outlined above, the odd-even 
staggering of experimental masses is not a pure measure 
of pairing correlations, but also has non-negligible contri- 
butions from the response of the underlying mean field 
to the blocking of a single-particle state. Here we can 
take advantage of the fact that various difference formulae 
are conceivable and take the recipe which best decouples 
mean-field and pairing properties. 

The odd-even staggering needs to be deduced from 
energy systematics. To that end, there are several finite- 
difference formulae in the literature to calculate A q from 
binding energies of adjacent nuclei. All available finite 



The number of the other kind of nucleons is assumed to 
be even and the same for all terms. Combining the ex- 
pansion (||^) of several adjacent nuclei leads to energy- 
difference formulae which can be used to approximate the 
gap Z\( b ). The two-point (first-order) formula leads to the 
one-nucleon separation energy, which mixes mean-field, 
single-particle and pairing effects strongly and should bet- 
ter not be used to fit the pairing strength. The next higher- 
order is the three point difference 



E(N + 1) 



ON 2 



2E(N ) + E(Nq 
1 <9 4 £ 



1) 



N 



12 ON 4 



(23) 



N 



+ ■■■ + D(N Q + 1) - 2D(N ) + D(N - 1) 



Assuming that the gap D varies only slowly with nucleon 
number and that the remaining contribution from the sec- 
ond derivative of Eq is negli gible (which is not so well ful- 
filled in some cases, see Sect. 4.2 and |52jj53|) this equation 



can be resolved into an approximative expression for the 
pairing gap 

A q 3) (N ) := ^ [e(N - 1) - 2E(N ) + E(N + 1^4) 

where 7Tjv = (— 1)^° is the number parity. A^ calculated 
from pure HF states without pairing but consideri ng p o- 
larization of the mean field is discussed in Refs. p2] , p3| in 
great detail. 

The next order corresponds to a four-point difference 
formula, but this order, employing an even number of nu- 
clei, gives an expression which is asymmetric around the 
nucleus with Nq and therefore offers two choices. With 
the same assumptions used going from (|2^) to (|I|), one 
possibility for the four-point gap is 



E{N - 2) - 3E{N - 1) 
+3E(No) - E{N + 1)1 . (25) 



This is an approximation for the gap at Nq — 1/2. The 
other possible four-point formula gives the gap at Nq + 1/2. 
The lowest-order derivative of Eq hidden in the four-point 
formula is now of third order. 

Although this four-point definition is widely used in 
the literature |8pq,p7|| we prefer the five-point formula 



difference formulae for A q are derived from the Taylor ex- Af\N ) := E(N + 2) - AE(N Q + 1) + 6E(N ) 



mansion of the nuclear mass in nucleon- 
56l 



number differences 



-4E(Nq - 1) + E{N - 2) 



(26) 



E(N) 



1 d n Eq 



i\ ON* 



(N - N ) n + D(Nq) (21) 



No 



where Eq is defined in ( fL9| ) and the Gap D given by 

even proton and neutron number, 
D = { A n odd neutron number, (22) 
zip odd proton number. 



which is symmetric and yields the best decoupling from 
mean field effects as we will see. The smooth contributions 
from the mean field to the gap are further suppressed com- 
pared to the lower-order formulae, the remaining deriva- 
tive of Eq entering A^ is of fourth order. 
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3.3 Quasiparticle Energies 

Another widely used approximation for the pairing gap is 
to calculate the energy difference (S3) by constructing the 
blocked many-body wave function of an odd-mass num- 
ber nucleus non-self-consistently from its BCS ground 
state as a so-called single-quasiparticle excitation . 
The lowest single-quasiparticle energy 



min(£' fc 



(27) 



- which we will simply denote as "quasiparticle energy" 
-E-quasi in the following - is then another approximation 
for the odd-even staggering. In the LN scheme the E). are 
given in first-order approximation by 



E k 



A 



2,8 



(28) 



see Appendix g for details. The important difference be- 
tween -Equasi and A^ is that for -E qua si the blocked many- 
body wave function is not calculated self-consistently. 



3.4 Spectral Gaps 

The calculation of Z\( 5 ) from mean-field models is a bit 
cumbersome since it requires information on five nuclei 
including nuclei with odd mass number. Moreover, the 
definition becomes inapplicable to describe the variation of 
pairing correlations with deformation for a given nucleus. 
At this point, we could recur to the purely theoretical 
definition ( p0| ) involving a blocked and an unblocked BCS 
calculation. This still requires involved calculations and 
can become unwieldy in deformed calculations. 

Therefore, a commonly used approach is to estimate 
the pairing gap from spectral properties of a nucleus. In 
schematic pairing models using the same pairing matrix 
element for all states the single-particle gaps ( |14| ) turn 
out to be state-independent and are thus immediately a 
measure for the pairing correlations. With local forces as 
used here we obtain state-dependent single-particle gaps 
Ak and have to define an average gap as representative 
for the strength of the pairing correlations. The authors of 
|l2| have proposed to use the average of the single-particle 
gaps (EP 



weighted with the occupation probability v\ 



(v 2 A) q 



(29) 



This definition, however, puts too much weight on deeply- 
bound states whereas pairing is a mechanism most active 
near the Fermi surface. We therefore propose an average 
with the same factor v^Uk as it appears in the accumu- 
lation of the pair density x( r )j see Eq. ([)]), yielding the 
spectral gap as 



(uvA) q = 



/ ken, 



fkV k u k A k 



J2kef2„ fk u kVk 



(30) 



Table 1. Pairing strengths V n for the neutrons and V p for the 
protons in MeV fm 3 for the DF and DDDI pairing functionals 
used in the BCS and LN schemes in connection with SkI4. 





DF 


DDDI 


Scheme 


V n Vp 


V n Vp 


BCS 


-323 -310 


-999 -1146 


LN 


-318 -250 


-947 -947 



Note that the spectral gaps are an estimate for the pair- 
ing gap and therefore the contribution of pairing correla- 
tions to the odd-even st agg ering (20). Assuming that the 
spectral gaps ( |2^ ) and (|30|) are an approximation for the 
square-root term in (|28| ) , approximately particle-number 
projected spectral gaps are given by 



(v-A)^ = {^A) q + \ 2s 
(uvA)f^ = (uvA) q + X 2 , q 



(31) 
(32) 



which was proposed by the authors of J46[ for the aver- 
age gap (v 2 A)( LN \ These spectral gaps will be discussed 
and compared with other alternatives for the calculated 
pairing gap in Section || 



4 Results and Discussion 
4.1 Fit of Pairing Strength 

The first step is to determine appropriate pairing strengths 
Vo }q for the DF and the DDDI functionals. We do that on 
the grounds of the five-point gap A^ and adjust the pair- 
ing strength by fitting calculated values for A^ to exper- 
imental ones for a large set of semi-magic nuclei, i.e. the 
isotope chains |fCa, ^Sn-^Sn, and gJPb-JJgPb, for 



the neutrons and the isotone chains IfCr 



i:s<> 



Xe 



_147 
54 82 



Th 



65 



and ?26 p °84 



215 
"126 



Ge 



94 



32-50R.U44, 



Acsg for the protons. A 



seperate fit has been performed for each one of the pairing 
functionals, DF or DDDI, and for each appraoch, BCS or 
LN. Each nucleon sort, proton or neutron, aquires its own 
pairing strength adjusted to isotonic chains, or isotopic 
chains respectively. 



The experimental data are taken from 1 58 . The result- 
ing values for the pairing strength are listed in Table [|. We 
obtain a reasonable fit of the pairing gaps for all pairing 
schemes, see also Figures |] (compare "Expt. with A^) 
and H in what follows. We have checked that the pairing 
strengths are not significantly changed (i.e. on the order of 
1%) when fitting instead theoretical A^> to experimental 
Z\ (3) or similarly the Z\ (4) . 



4.2 Comparison of the Finite-Difference Formulae 

The various finite-difference gaps were introduced as ex- 
perimentally accessible approximations to A^°\ i.e. the 
energy differences of blocked and unblocked calculations 
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Fig. 1. Comparison of the three-point A^ 3 \ four-point A^ and five-point difference formulae A^ and the blocking gap 
A (h \ The left panel shows results for neutron gaps in the chain of tin isotopes, the right panel for proton gaps in the chain of 
the N — 82 isotones. 



of odd-mass nuclei. Figure [I] compares directly the perfor- 
mance of the Z\W, i = 3,4,5, in this respect for calcula- 
tions within the BCS approach. The three-point gaps A^ 
show a pronounced odd-even staggering for all N = 82 
isotones and the tin isotopes with neutron numbers below 
the TV = 82 shell closure. Note that we have disentangled 
that by drawing a separate line for even-even and odd- 
mass nuclei. These two lines for A^ embrace the gap 
Z\( b ). The stag gering is caused by non-vanishing mean- 
field contributions to Z\^ 3 ^ (mainly the second derivative 
term in Eq. (|23|)), which enter with a different sign for 
even-even and odd-mass nuclei. The four-point gaps A^ 
from Eq. ( |25| ) give a smoother approximation for A^ h \ 
But as expected the A^ are slightly shifted versus the 
/\( b ) which becomes rather obvious where the gaps change 
rapidly, e.g. around N = 63, N = 70, Z = 57, Z = 65. The 
five-point gaps A^ give the best overall agreement with 
the A<- h \ 

The oscillation of the A^ around the values for A^ 
has a simple geometrical reason, as can be seen from Fig. |[ 
is per definition (|2C| ) the shift between the smooth 
curve connecting the (unblocked) BCS ground-state en- 
ergies Eq of all nuclei (thick line through circles) and 
the smooth curve that connects the (blocked) ground- 
state energies E\> of odd-mass nuclei (dotted line through 
full boxes). The three-point gap of an even-even nucleus 

4vL(iVo) = ±[E h {N - 1) + E h (N + 1)] - E (N ) is the 
shift between the smooth curve connecting the Eq and 
the average energy of the two adjacent odd-mass nuclei, 
which lies outside the band given by Z\( b ). From this fol- 
lows immediately Zicvcn > A^ for bound nuclei. A simi- 
lar construction leads to A^J d < A^ h \ see Fig. |J For Z\( 4 ' 
and Z\( 5 ) the deviation from A^ becomes of course much 
smaller because the higher-order finite-difference formu- 
lae give a better approximation of the smooth curves con- 
necting the Eq and Et, respectively. This qualitative re- 
sult does not depend on the level of sophistication for the 
calculation of the odd-mass-number nuclei. Considering 



polarization effects will shift E^ with respect to Eq (which 
has to be counterweighted by a refit of the pairing strength 
[ p4|j55|] ) and may distort the surface of the E\>, but will not 
change the sign of its curvature. 

There remains a significant difference between all A^ l \ 
i = 3, 4, 5 and A^ around shell closures where finite- 
difference formulae (except A^ for odd nuclei) produce a 
peak which becomes broader with increasing order of the 
difference formula. We want to discuss the origin of this 
peak for the example of the five-point gap A^ . The five- 
point approximation d26| ) for A^ makes three assump- 
tions: (i) the binding energy is an analytical function of 
the nucleon numbers and therefore can be expanded in a 
Taylor series in the range of two mass units around the 
considered nucleus; (ii) derivatives of Eq of higher than 
third order are negligible; and (iii) the gap varies only 
slowly with nucleon number. The first two assumptions are 




Particle Number 



Fig. 2. Schematic comparison of A^ and A^ . Eq and E b de- 
note BCS ground-state energies and the energy of the blocked 
ground state of an odd-mass nucleus respectively. Filled sym- 
bols denote binding energies that are experimentally accessi- 
ble, while the BCS ground state of an odd-mass nucleus (white 
circles) can be calculated from the mean-field model only. 
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strongly violated at shell closures, where the kink in the 
systematics of binding energies does not allow the Taylor 
expansion (|2l]) and leads to a spurious contribution from 
the mean-field functional to the finite-difference gaps. 
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formulae like A^ ss 12/v^4 are intended to represent the 
average trend of the odd-even staggering, not the pairing 
gap. The discrepancies between the two quantities can be 
expected to decrease with increasing mass number and to 
be very small for heavy nuclei. 

An analysis of gaps from a complementary point of 
view is given in ]52j for the case of A^ 3 \ This study omits 
pairing altogether and concentrates exclusively on polar- 
ization effects for A^ using pure deformed HF calcu- 
lations. The focus is on small nuclei because these have 
most pronounced deformation effects. In this framework, 
it is found that 3-point gaps show a pronounced odd- 

(3) 

even staggering with the ^4 odd being close to zero while 

(3) 

most of the Z\cvon have large positive values in most cases. 
This finding is explained in terms of the macroscopic- 
microscopic model in that a large contribution from the 
symmetry energy is counterweighted by the difference of 
single-particle energies when calculating A oc ( d . This ob 



Fig. 3. Five-point gaps Zr 5 ' computed using either blocked 
calculations of the odd-mass nuclei (dashed lines) or their BCS 
ground states (dotted lines) for the chain of tin isotopes with 
BCS+DF pairing. Non-zero Z\' 5 ' from not blocked calcula- 
tions appear only at closed shells and are caused by a spurious 
contribution of the mean field to /A (5) . The difference of both 

curves (solid line) gives the contribution from the blocking to 
4 (5)_ 



To visualize this effect wc compare in Fig. | A^ with 
five-point gaps Z\( 5 )( nb ) using the binding energies of not 
blocked calculations of the odd-mass-number nuclei, which 
carries only the the spurious mean-field contribution to 
A^. Subtracting Z\( 5 )( nb ) from A^ one gets the contri- 
bution from the blocking of the odd particle to Z\' 5 ' . The 
peaks at closed shells disappear yielding a smooth curve 
for Z\( 5 ) - Z\( 5 )( nb ) which now follows the values of A^ 
everywhere. 

Besides the immediate vicinity of shell closures, wc find 
that A^ is a very good approximation for the staggering 
gap Z\( b ) . For a few additional nuclei there remains a small 
difference between A^ and A^ b \ see Fig. |3|. This occurs 
for example around N — 68 and N — 72 for the tin iso- 
topes. For those nuclei the assumption of a slow variation 
of the gap with nuclcon number which enters the five- 
point formula (26) is not valid: the gap changes in these 
regions by about 40% (due to a sudden change of the den- 
sity of single-particle levels at the Fermi surface in these 
nuclei). But even then the deviation between A^ and 
/\(b) remams acceptably small. We thus prefer A^ for 
the fit of pairing strengths. The small difference between 
Z\( 5 ) and Z\( b ) even suggests a simplified fitting procedure 
where calculated A^ 1 are adjusted to experimental values 
for Z\( 5 ) in chains of spherical semi magic nuclei. 

As explained in Sect. 3.1, the A^ unavoidably con- 



tain a contribution from the polarization of the mean field 
in odd nuclei. It is therefore somewhat unlucky and con- 
fusing that Z\( b ) is usually denoted as "pairing gap" in 
the literature, see e.g. [pl| , p0|j57|] . The commonly used fist 



servation led the authors of 
v(3) 



to the conclusions that 



Z\ odd is very close to the pure pairing gap while Zieven 
contains a contribution from the mean field. The higher- 
order gaps Z\( 4 ) and A^> turn out to be rather useless in 
that environment. We want to point out that the odd-even 
staggering of Z\( 3 ) observed in [[52) is not the staggering of 
Z\( 3 ) around the values for A^> as we discuss it here, see 
e.g. Fig. [|. The deformation staggering is a phenomenon 
much similar to the spurious peak of finite-difference gaps 
at major shell closures discussed above. Mind that pure 
deformed HF calculations in small nuclei produce a sub- 
shell closure for each even-even nucleus by virtue of the 
Jahn-Teller effect. This, in turn, leads to a spurious con- 

(3) 

tribution to zicven for nearly all even-even nuclei which is 
half of the energy difference between single-particle levels 
[ p2[ , and it has devastating consequences for the system- 
atics of 4-point and 5-point gaps. The picture changes 
dramatically when pairing is included, as we do here. Pair- 
ing induces a drive to spherical shapes and thus reduces 
deformation effects dramatically while rendering blocking 
the the dominant contribution to the gaps. This smoothes 
the systematics of binding energies and thus of any 
Moreover, we are discussing here semi-magic medium and 
heavy nuclei which have spherical BCS ground states. We 
are thus considering a sample with minimal mean-field ef- 
fects, just appropriate to concentrate on pairing features. 
The interplay of polarization effects which we neglect here 
and pairing correlations will probably play a role for small 
non-magic nuclei. It deserves further inspection in the fu- 
ture. 



4.3 Spectral Gaps in BCS Pairing 

Having discussed the various finite difference gaps, we now 
take the five-point gap as reference value and study the re- 
lation to the spectral gaps and quasiparticle energies. Fig- 
ure |^ compares calculated results for A^ 5 \ (v 2 A), (uvA) 
and -Equasi computed within the BCS approach. 

Let us start the discussion by looking at the neutron 
gaps in the chain of tin isotopes calculated with the DF 
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Neutron Number N Proton Number Z 

Fig. 4. Comparison of the five-point gap A^\ the spectral gap (uvA), the average gap (v 2 A), and the single-quasiparticle 
energy -Bquasi with experimental values for the five-point gap for tin isotopes (left panels) and N = 82 isotones (right panels) 
calculated with a delta pairing force (upper panels') and DDDI force (lower panels) in the BCS scheme. The experimental values 
are calculated from the binding energies given in jig] using the five-point formula A^ (^fi|). 



pairing interaction in the BCS approach (upper left panel) . 
The spectral gaps (uvA) and (v 2 A) show a pronounced 
odd-even staggering where the gaps of even-even nuclei 
have larger values than those of odd-mass nuclei. This is 
caused by the blocking of one state with a large weight uv 
in the odd-mass nucleus. The blocked state does not con- 
tribute to the pairin g po tential, leading to overall smaller 
single-particle gaps The amplitude of the odd-even 
staggering decreases, of course, with increasing neutron 
number TV because the relative contribution of a particular 
state to the pair density becomes smaller with increasing 
level density in the heavier isotopes. 

Both (v 2 A) and (uvA) are exactly zero for closed shell 
nuclei, i.e. N — 50, N = 82, N = 112, and the adjacent 
odd mass-number nuclei. In these nuclei the BCS scheme 
breaks down. This is a deficiency of the BCS scheme which 
is related to the particle-number uncertainty of the BCS 
state @. 

The quasiparticle energies -©quasi show a similar depen- 
dence on the neutron number as the A^ , but for most 
nuclei they are 100-250 keV larger (and thus the same 
amount larger with respect to the A^). This is caused by 
calculating the blocked many-body wave function entering 
©quasi not self-consistently. The variational principle be- 
hind the self-consistent calculation of the odd-mass nuclei 



entering the A^ leads always to larger a binding energy 
of the odd-mass nuclei, lowering the calculated gap. 

The quasiparticle energies show the same peak at shell 
closures as the finite-difference gaps which is again related 
to a spurious contribution from the mean field. For non- 
magic nuclei the lowest single-quasiparticle state usually 
corresponds to a single-particle state with e' k f=a X q leading 
to the single-quasiparticle energy -E qU asi ~ + ^2,q- In 
magic nuclei one has an additional contribution from the 
first term in the square root in Eq. ( p8| ) since the Fermi 
energy is approximately in the middle of the gap in the 
single-particle spectrum (In the BCS scheme, where the 
pairing breaks down for closed-shell nuclei the derivation 
of Eq. ( p8| ) is not valid. Then one has different Fermi ener- 
gies for the removal and addition of a particle, which are 
the single-particle energies of the last occupied and first 
unoccupied single-particle state respectively). The large 
jump in the Fermi energy is the reason for the kink in the 
systematics of binding energies at shell closures, which in 
turn causes the peak in the finite-difference gaps discussed 
above. 

While all definitions of the gap give similar values in 
the valley of stability, there appear large differences be- 
tween the A^ and -E qua si on one hand and the (v 2 A) and 
(uvA) on the other hand for neutron-rich nuclei beyond 
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the N = 82 shell closure. The spectral gaps overestimate 
the A^ (which closely follow the A^ as discussed ear- 
lier). The (uvA) are smaller than the (v 2 A) for all tin 
isotopes and in most systems the (uvA) are closer to the 
Z\< 5 > than the (v 2 A). 

Now we want to look at the changes when employing 
the DDDI pairing functional instead of the simpler delta 
force, see the lower left panel of Fig. ||. There are two sig- 
nificant differences to the results obtained with the DF 
functional: (i) all gaps except (v 2 A) are larger by roughly 
20 % for neutron-rich nuclei with N > 82, and (ii) the 
(v 2 A) are always much smaller than the (uvA) . This trend 
is just the opposite from the one for the DF functional. 
This is caused by the different choice of weights in the def- 
inition of (v 2 A) and (uvA) in combination with the spec- 
tral distribution of the single-particle gaps Ak- While the 
DF pairing potential follows roughly the nuclear density 
distribution, the DDDI functional gives a pairing poten- 
tial which is sharply peaked at the nuclear surface, see 
Fig. ||. This leads in case of the DF interaction to single- 
particle gaps Ak of comparable size for all bound states 
while in case of the DDDI interaction the Ak of deeply 
bound single-particle states are rather small, see the mid- 
dle panel of Fig. ^. Together with the weight factors used 
to calculate (uvA) and (v 2 A) - see the upper panel in 
Fig. ^ - this gives the observed pattern for the spectral 
gaps: when calculated with the DF interaction they are 
rather insensitive to the choice for the weight factors while 
there is a huge difference in case of the DDDI interaction. 

This explains also why the A^ extrapolate quite dif- 
ferently when comparing DF and DDDI pairing for neutron- 
rich nuclei. In these nuclei the states at the Fermi surface 
are only loosely bound and therefore have a large spatial 
extension but only small overlap with the volume-like DF 
pairing potential. An extreme example is the tin isotope 
with N = 112 where DF pairing breaks down but DDDI 
pairing is still fully active. Experimental data on the exci- 
tation spectra of neutron-rich nuclei will give in the future 



valuable information to distinguish between volume-like 
(DF) and surface-peaked pairing interactions. 
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Fig. 6. Weights v 2 (full diamonds) and 2uv (open diamonds) 
entering the calculation of the spectral gaps (v 2 A) and (uvA) 
(upper panel), single-particle gap Ak in MeV (middle panel) 
and single-particle root-mean-square radius (lower panels) of 
the neutrons in 124 Sn, calculated with SkI4 in the BCS scheme. 
Weights and radii are shown for results with the DF functional 
only (DDDI gives very similar results) , while the Ak are shown 
for the DF (diamonds) and DDDI (circles) pairing functionals. 
The horizontal line in the lower panel is the root-mean-square 
radius of the neutron density. The uv in the upper panel are 
multiplied by two to have the same possible maximum value 
of one as the v 2 . 
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Fig. 5. Local pairing potential A(r) of the neutrons in 124 Sn 
for the DF and DDDI pairing functionals. The DF pairing po- 
tential acts over the whole volume of the nucleus, while the 
DDDI pairing potential is peaked at the nuclear surface. 



This different behavior of (v 2 A) on one hand and (uvA) 
and Z\( 5 ) on the other hand hints that is possibly danger- 
ous to use different definitions of the gap for experimental 
and calculated gaps when fitting the pairing strength and 
comparing calculated and experimental values. 

The right panels of Fig. [| show the pairing gaps of the 
protons in the chain of N = 82 isotones. Qualitatively the 
results are similar to those for the neutron gaps in the tin 
isotopes. But the overall reproduction of the experimental 
values is much better than in the case of tin isotopes and 
the differences between the various gaps is much smaller 
when going towards the drip-line. The Coulomb poten- 
tial stabilizes even loosely-bound protons. Therefore the 
difference between the gaps comparing volume-like and 
surface-like pairing potentials is quite small. Only for the 
DDDI force remains the large difference between the spec- 
tral gaps (v 2 A) and (uvA) explained above. 
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4.4 Lipkin-Nogami Pairing 

Figure shows the gaps of the neutrons in the chain of tin 
isotopes, now calculated with the LN scheme. The spec- 
tral gaps (v 2 A) are omitted here. Instead we compare the 
"bare" (uvA) (|30[) with the particle-number corrected val- 
ues (uvA) + A2 (|32|). 

The global pattern of the various gaps looks very sim- 
ilar to the one obtained with the BCS scheme, see Fig. [|. 
The most obvious difference between the BCS and LN 
methods is that the LN scheme does not break down 
for closed-shell nuclei. Therefore the "bare" spectral gap 
(uvA) has a finite - but still somewhat too small - value 
around magic nuclei. Adding A2 gives better results around 
shell closures, but away from shell closures the difference 
between (uvA) and (uvA) + A2 for the tin isotopes is too 
small to decide on one preferred definition. 

From the difference between (uvA) and (uvA) + A2 it 
can be seen that A2 is largest around shell closures. But 
this indicates also that the LN approximation might not 
be sufficient for magic nuclei, a variational calculation of 
A2 or even full projection of the many-body wave function 
is needed |2(J there. 

The single-quasiparticle energies i? qua si follow closely 
the particle-number projected (uvA) p + \2, P , which is eas- 
ily understood remembering that A2 is added to both 
quantities. At shell closures, however, the single-quasipar- 
ticle energies overestimate the experimental gaps. Like in 
the case of the BCS scheme the -E quas i are nearly always 
larger than the calculated five-point gaps A^ 5 \ 

4.5 Comparison of all Models 

Finally, we want to compare the four pairing models, i.e. 
any combination of BCS or LN and DF or DDDI pairing. 
The comparison is done with respect to the five-point gap 
Z\( 5 ), which we prefer as the most robust empirical defini- 
tion and which has turned out to be the most useful defi- 
nition of the pairing gap. In Fig. ^ the five-point gaps cal- 
culated from various pairing schemes are compared with 
experimental values. The differences between BCS and LN 
are generally very small. It is to be remembered, how- 
ever, that the effective pairing strength is readjusted for 
the LN scheme. The (approximate) particle-number pro- 
jection increases the total binding energy, but this effect 
is renormalized by virtue of the fit delivering a slightly 
smaller pairing strength in the LN scheme. There is one 
detail where BCS and LN differ: the peak of A^ in the 
vicinity of closed shells is more spread out in case of LN 
which is probably due to the softening of the shell closure 
byLN. 

The differences between the pairing forces (DF versus 
DDDI) are much larger. For the neutron gaps of tin iso- 
topes close to the valley of stability and the proton gaps 
in the N = 82 isotones all schemes and forces still give 
similar results, but large differences between the DF and 
DDDI force occur for neutron gaps around N = 60 and 
very neutron-rich nuclei with N > 82. Only the DDDI 
model can describe the gaps in both the light and heavy 



known tin isotopes, while the DF interaction overestimates 
the gaps in the light ones by up to 15-20%. The particle- 
number projection has only a small effect on the gaps when 
the strength of the pairing interaction is readjusted. The 
better description of the tin isotopes around TV = 60 with 
the DDDI force gives a hint that this pairing interaction 
may be more realistic than a delta force. This, however, 
has to be taken with a grain of salt: it may also be a spuri- 
ous effect due to a deficiency of the underlying mean-field. 
The disagreement between calculated and experimental 
Z\( 5 ) around N = 70 is rather robust in that all forces and 
schemes give very similar results. For a profound decision 
which pairing interaction gives the most realistic results 
throughout the chart of nuclei more and other observables 
have to be investigated using various forces for the under- 
lying mean-field. Research in that direction is underway. 

We have checked that one obtains very similar results 
when comparing experimental and calculated three-point 
and four-point gaps. However, it is important that the 
experimental and calculated values are computed from the 
same formula. As already shown in Fig. |l|, the various 
finite-difference formulae may give results which differ by 
25%. 



5 Summary 

We have compared various approximations for the cal- 
culation of the pairing gap: three-point, four-point and 
five-point finite difference formulae where the gap is esti- 
mated from total binding energies of adjacent nuclei, spec- 
tral gaps with different weights which put bias on well- 
bound states or levels at the Fermi surface, and the single- 
quasiparticle energy. Predictions of four different pairing 
models for the gaps were compared, namely the BCS and 
LN pairing schemes employing the DF or DDDI interac- 
tions. Experimental values for the pairing gaps are usually 
calculated from a finite-difference formula. For the calcu- 
lated gaps we find that apart from shell closures there 
are non-negligible deviations up to 25% between the vari- 
ous definitions. Some definitions of the pairing gap cannot 
be used for closed-shell nuclei. Therefore, it is the safest 
choice to compute the pairing gaps from mean-field cal- 
culations in the same way as the experimental values. 

The natural definition for the calculated gap is the 
difference in binding energy A^ h>> between the fully paired 
BCS ground state and the blocked one of odd-mass nu- 
clei. The five-point gap A^ provides a reliable quantity 
to fit the effective pairing strength, among all approxima- 
tions for the pairing gap it is closest to A^ h \ Therefore we 
use it to calculate the experimental pairing gaps and take 
it as point of reference for all other approximations for 
the pairing gap. Like the other finite-difference formulae 
contains a spurious contribution from the mean field 
around magic numbers, which is related to the jump of 
the Fermi energy at shell closures. 

The four-point gap A^ has nearly the same overall 
quality as the five-point gap as compared to A^ h \ but its 
definition has an ambiguity, so that we prefer the five- 
point gap. The three-point gap A^' has a large contribu- 
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Neutron Number N Proton Number Z 

Fig. 7. Comparison of the five-point gap A^\ the spectral gap (uvA), the particle-number corrected spectral gap (uvA) + A2 
and the single-quasiparticle energy -E qU asi with experimental values for the five-point gap for tin isotopes (left panels) and 
N = 82 isotones (right panels) calculated with a delta pairing force (upper panels) and DDDI force (lower panels) in the LN 
scheme. 




tion from the mean-field, which becomes rather obvious 
looking at proton gaps, consequently this quantity should 
not be used for the fit of the pairing strength since deficien- 
cies of the underlying mean-field (especially concerning 
the symmetry energy) may be visible in the A^ . When 
comparing calculated and experimental A^ respectively 



Z\( 4 ) for a well-adjusted interaction, however, both gaps 
show the same quality like the five-point gaps A^ 5 \ 

In situations where it is not possible to calculate A^ , 
e.g. looking at potential energy surfaces, a reasonable ap- 
proximation for the pairing gap is provided by (uvA) + A2 . 
By adding A2 to the spectral gap the approximate particle- 
number projection gives an improved reproduction of the 
calculated A^ in most cases. The spectral gap (uvA) + X 2 
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is in much better agreement with the calculated A^ than 
the averaged gap (v 2 A) + A2, which sets too much bias on 
deeply bound states and therefore should not be used in 
models in which the pairing interaction acts mainly at the 
nuclear surface. The single-quasiparticle energies always 
overestimate the A^ b \ but this is not too surprising be- 
cause they are calculated from a non-self-consistent wave 
function of the odd nucleus. 

The results show the danger of comparing calculated 
and experimental gaps computed from different defini- 
tions. The deviations between the various definitions de- 
pend on the actual nucleus and become generally larger 
when going towards the drip-lines. This is important es- 
pecially when adjusting the free parameters of the pairing 
interaction. Experimental and calculated gaps should be 
calculated in the same manner, the five-point gap provides 
a useful tool for that. 

The choice of the test cases (heavy semi-magic nu- 
clei) and restriction to spherical nuclei have minimized the 
impact of polarization effects on the gaps. Recent explo- 
rations of dynamical polarization effects 1 54 5f| hint that 
these may be non-negligible at a quantitative level. Con- 
clusive answers are yet inhibited due to uncertainties of 
present mean-field models in the time-odd channel. This 
point deserves attention in future work. 

Another open question is the functional form of the 
pairing interaction. We find in this paper that the DDDI 
functional gives a slightly better description of pairing 
gaps compared with a delta pairing force, but to give a 
definitive answer one has to look at more nuclei and other 
data as well. Work in this direction is underway. 
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A The calculation of A 2 
A.l The General Expression 

In case of an underlying many-body Hamiltonian H the 
parameter A2 in the variational equation is fixed by the 



condition 

where K is given by 

k = h - 



(KN. 



]T (\i, q (N q ) + \ 2 , q {N q 2 ) 

q£p,n 



(33) 



(34) 



and N2 is the two-quasiparticle part of the particle-number 
operator ( |l7j ) . Because we look at like- particle pairing only 
the equations for protons and neutrons separate. There- 
fore the index q for the isospin of the particle-number op- 
erator, the single-particle states etc can suppressed to get 
a compact notation. Introducing the shifted many-body 
state |B8l 



10 



|0) 



with 



|o> = 10 



£=0 



the condition (|33|) is equivalent to 

df<o|£|o| s=o = o . 

This can be resolved into an expression for A2 
d 2 A0\H\O 



(35) 



(36) 



(37) 



£=0 



The denominator of 
Wick's theorem 



n) can be calculated easily using 



a £ 2 (o|iv 2 |0 



£=0 



fc3=o 



2 2 
U k V k 



16 J2 u i v i -( 38 ) 



So far we have formulated the numerator in terms of an 
underlying many-particle Hamiltonian, but the formula- 
tion (|36| ) has the advantage that it can be translated into 
the formal framework of energy functionals p8j 



<o|#IO 



£>(€) 0*(0l 



(39) 



p(® and are the shifted density matrix and pair den- 
sity matrix respectively, which have to be calculated now 
as non-diagonal matrix elements 









= x (0( 


r (£) 





(x,x') 
(x,x') 



)|^ t (x')V'(x)|0 (40a) 
-<7'(0|^(r',-</>Kx)|0 (40b) 
-<j(0\ft(x')ft(r,-a)\O (40c) 



All local densities and currents the Skyrme and pairing 
energy functionals depend on can be derived from these 
density matrices. The second derivative of the energy func- 
tional is given by 

-trtH - - d 6 pi d i p 2 + 2 J .^ it d£XO£X 



SpiSp 2 



Sx^x* 



5 2 £ 8 2 £ 1 

-2— d i pda + 2——d t pd i r ■ (4i) 



SpSx 



SpSx* 
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6 stands for a functional derivative. The trace is a short- 
hand notation for the integration and summation over all 
coordinates 



(68 . 



dx dx' 



5£ 



5p(x,x') 



7TP( X ' X ') 



(42) 



The terms with a single trace in (|4lj) vanish in case of a 
BCS ground state, while the terms with double traces can 
be simplified defining the response density matrices p, x 
and x* 



p = — i p^ 



X = -i X (4) 



(43) 



leading to the final expression 



d 2 £^ 



-tr tr 



'Pi P2 



2 — — , X X 



<^X* 

n <5 2 £ . . „ 5 2 £ _ _„ 
+2 T ^rPX + 2 7 ^ T — PX 



(5/5 (5x* 



(44) 



which has to be inserted into (|37|). Since we are consider- 
ing pairing between like particles only there are no mixed 
terms with derivatives with respect to proton and neu- 
tron densities. This expression has to be evaluated now 
for the pairing and the mean-field energy functional. This 
density-functional approach to the calculation of A2 in- 
corporates in a natural way the additional contributions 
to the LN equations arising for density-dependent inter- 
actions discussed in [|47|j4S 



A. 2 The Linear Response of a Skyrme Energy 
Functional 

The Skyrme energy functionals are constructed to be effec- 
tive interactions for nuclear mean-field calculations. For 
even-even nuclei, the Skyrme energy functional used in 
this paper 



£ — £kin + £sk + £1 



(45) 



is the sum of the functional of the kinetic energy £kin, the 
effective functional for the strong interaction £sk and the 
Coulomb interaction £q including the exchange term in 
Slater approximation. The actual functionals are given by 



£kin 

£c 
£sk 



2m 



d A r t 



d A r 



d 3 r d 3 r 



2 H 



3 /Pp( r )Pp( r ') 



r — r 



3e* 
4 



1/3 



b\pr pAp 



+2 



E 



b' 



2 

b'lPqTg 



d 3 rpi /3 



b A pV ■ J 



b' 



'-PqApq 



(46) 



The local density p q , kinetic density T q and spin-orbit cur- 
rent Jg entering the functional are given by 



Pq(*) 



Tq(r) 



ken, 



p g (r,fj;r,cr) 
V-V'p g (r,(r;r» 

£«*|V0 fc (r)| a , 

i(V-V') x Y Pq(r,a;r',a')a a , a 

a, a' — ± 

1 £ 4 kWVxa^(r)- h.c 



keO„ 



(47) 



with q £ {p, n}. a a i a is the matrix element of the vector 
of the Pauli spin matrices between the unit spinors with 
spin projection tr'/2 and a/2. Densities without index in 
( ft6| ) denote total densities, e.g. p — p p + p n . The (f>k are 
the spinors of the single-particle wave functions, the v 2 
occupation probabilities. The parameters bi and b' i used 
in the above definition are chosen to give a most compact 
formulation of the energy functional, the corresponding 
mean-field Hamiltonian and residual interaction [ ]59[ . The 
Skyrme energy functional contains an extended spin-orbit 
coupling with an explicit isovector degree-of-freedom as 
used in the parameterization SkI4 p4j . 

The response densities needed for the evaluation of 
(Eih are given by 



2 J2 44 \M 



k£f2„ 



fq=2 E 44 |V0fc| 



ken„ 



J 9 = -i E 4 



h\. V X &(j) k - (V X a0 fc )Vfc 



ken„ 



Evaluating Eq. (Q) for the Skyrme energy functional 
leads to 



€=0 
tr tr|pg : i 

+2p q 



S 2 £ ; 



Sk 



Sp q .l 5p q ,2 

S 2 £ Sk 

Sp q SJ q 



Pq.2 + 2p q 



S 2 £ ; 



Sk .. 



Sp q 5r q q 



d 3 r {(b Q -b' )p 2 q + 2(b 1 -b' 1 )p q fq 

~{b 2 - b' 2 )p q Ap q - 2(64 + b' 4 )p q V ■ J q 
+ i[(a + 2)(a + l)b 3 -2b' 3 ] p a p\ 
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For the local densities appearing in ( fig ) the trace reduces the local density 
to a spatial integral. For the protons one has an additional 
contribution from the Coulomb interaction 



n2 e«) 



tr tr <^ p p 



8 2 £ c 



Pp 



d 3 r p p V cou i 



1/3 



e 



3 
n 

Z 2 



1/3 



d 3 r Pp 



2/3 
Pp 



Z 2 



d 3 r Pp 



2/3 
Pp 



(49) 



Poisson's equation for the response Coulomb potential 



, , in- -' C-DDDI ~* , „_ CDDDI - 

_ <5 2 £dddi - 



<v 9 <>xS q 



T /d 3 r<jx* X, 



<5 2 £ 
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Pq,2 
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(54) 



AV C( 



-4tt p p 



(50) B Single-Quasiparticle Energies 



is solved numerically using the techniques explained in 
pp|l . The kinetic energy gives no contribution to A2. We 
omit the contribution from the center-of-mass correction 
(and therefore the approximate particle-number correc- 
tion of this term). 



A. 3 The Linear Response of the Pairing Energy 
Functional 



For the calculation of the contribution of a pairing energy 
functional of type (|J) to Eq. ( f44| ) the local response pair 
density \ q is needed 



X q = -i d e xi° 



-4 J2 u k vk w 



(51) 



feen g 

k>0 



The response pair density is not hermitian, the adjoint 
response pair density x* reads 



Xq = "I 0$X q 



-i dtv*M 



4 fk u * v k 10*1 



(52) 



ten, 



In case of the DF pairing energy functional there is only 
one contribution from the derivatives with respect to the 
pair density 



{ DF 



-DF 



2trtr 1X 9 J^ZZ Xq 



$Xq $X*q 



"' 1 d 3 r x* X q 



(53) 



The single-quasiparticle energies are deduced from the 
non-self-consistent ansatz 



\k) 



&\ |BCS) 



v k a 



|0) (55) 



a>0 



for the blocked many-body wave function of the odd mass 
number nucleus |[|,^]] ■ The single-particle wave functions 
of all states and the occupation probabilities of all un- 
blocked states m k are taken from the self-consistent 
calculation of the BCS ground state. The normalization of 
|fc) requires v k = 1, u k = v k = v% = 0. Note that \k) does 
not yield the proper particle number of the excited state 
because the occupation of the blocked pair of states (fc, k) 
in the BCS ground state in general will differ from one. 

The excitation energy E k of the lowest one-quasipar- 
ticle state is an approximation for the odd-even mass dif- 
ference. To take the difference in particle number between 
the fully paired BCS state and |A;) into account, E k has 
to be calculated from the difference of K. given by ( |Tc| ) 
calculated for and the BCS ground state. This leads 
to 

E k = E {k) + (k\N'\k) - £<® - (N') (56) 
where we have introduced the abbreviation 



<j£{p,n} 



(57) 



£(°) is the energy functional of the BCS ground state, 
while £( fc ) is the energy functional evaluated for the one- 
quasiparticle state \k). 

The expectation values of N' are conveniently calcu- 
lated with standard operator techniques [ 



(58) 



(k\N'\k) = {a k N' &t) = <, + (a.N^al) 



In case of the DDDI pairing energy functional there are where Nq is the BCS expectation value of N' while 
additional contributions from derivatives with respect to is its (11) component. The Bogolyubov transformation of 
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the particle-number operator reads 



fcS=0 



[ v k + ~ v k) & l & k + u k v k {al&\ + a- k a k ) 

fc3=0 



= ^00 + ^11+^20 + ^1 



02 



(59) 



(N 2 )u can be calculated from ( p9|) without repetition of 
the quasiparticle transformation. One obtains 



(60) 



The calculation of the contribution of a many-body Hamil- 
tonian operator to E k can be found in many textbooks, 
see e.g. [|,^,||. Here, however, we want to calculate it in 
the framework of effective energy functionals. The energy 
functional = £[p {k) , x (fe) , X* (fc) ] of the one-quasiparticle 
state is to be calculated from the density matrix and pair 
density matrix evaluated for the one-quasiparticle state 
|fc) 

p( fc )(x;x')Hfc|^( x 'Wx)|fc) 
= ^(x')^(x)+^^^(x')0 m (x) , (61a) 



n>0 



,(*0 



; )(r) = -]>>(fc|V3(r,-a)V;(x)|fc) 



mVra 



(61b) 



n>0 



p( k \x;x') and X ( r ) are inserted into the definition of 
the energy functionals ( p5| ) and (0). This leads to 

E k = £( fe > - - (A, - 4A 2 ,,t£)(«2 - t£) + A 2 , g (62) 

where A g = Ai l9 + 4A 2i g(A r <j + 1) was used. 

The usual approximation for the one-quasiparticle en- 
ergy is obtained taking only the linear change in the den- 
sities into account 



6p 



15S_ 
2^ 



1 58 

25r 



with 



5p = pM-p={l-2vl) <f>Ux') 0fc(x) 
5 X = 5x* =X {k) -X = 2u k v k \M 2 ■ 
This yields 

£W _ £(0) ~ efc ( u | _ w 2) + 2/ fc Z\ feUfeUfe 



(63) 

(64a) 
(64b) 

(65) 



Together with the contributions from the particle-number 
operators, Eqns. (59) and (|60|), one obtains 



E k = (e' k - \ q ) (u% - vl) + 2f k A k u k v k + A 2 



(4-A g ) 2 + / fe 2 ^ + A 2>9 



(66) 



where the definitions of the Fermi energy (|T^), the renor- 
malized single-particle energy ( |l3| ) and the expressions for 
v\ and = 1 — v\ have been inserted. 
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